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We explore a novel method of mass reconstruction in events with missing transverse momentum 
at hadron colliders. In events with sizeable boost factors in the final steps of dual multi-stage decay 
chains, the missing energy particles may each be approximately collinear with a visible standard 
model particle, spanning a narrow "MET-cone." We exploit this collinear approximation, when 
applicable, to reconstruct the masses of exotica. 



Introduction — The start of the Large Hadron Col- 
lider (LHC) at CERN gives hope to the discovery of long- 
anticipated TeV scale new physics. In many of these 
scenarios, discrete symmetries are motivated by the re- 
quirement that higher dimensional operators which con- 
tribute to processes such as weak gauge boson form fac- 
tors and baryon number violation must be highly sup- 
pressed. Such discrete symmetries may simultaneously 
be responsible for stabilizing the dark matter compo- 
nent of our universe. The typical collider signatures of 
such discrete symmetries are characterized by events with 
large amounts of missing transverse momentum. The 
undetected particles in such events, the lightest parti- 
cle charged under the discrete symmetry, complicate the 
reconstruction of the masses of new exotica. This is par- 
ticularly the case at hadron colliders, where the initial 
parton momenta are unknown. There has been recent 
substantial progress in mass measurement in such sce- 
narios over the past few years (see [T] for a recent review, 
and complete citation list). Most of these methods re- 
lies on the kinematics of the events and fall into several 
broad categories: invariant mass endpoint methods [2H1], 
mass relation/polynomial methods [IHH], M^-like meth- 
ods [SHU], and various combinations of them. The ex- 
istence of multiple methods is crucial, providing comple- 
mentary techniques for extracting information about the 
underlying physics model. 

In this letter, we explore a conceptually new method 
of mass determination that is particularly useful when a 
decay chain terminates with the disintegration of a rela- 
tively boosted exotic particle to the lightest exotic plus 
a visible standard model (SM) particle. 1 Such events 
are characteristic of models in which pair produced color 
charged exotica are quite heavy, and in which there is 
simultaneously a small amount of phase space available 
for the NLSP decay. As a motivator for such scenarios, in 
the MSSM, the LEP II bound on the higgs mass prefers 



a large top squark mass, or in general TeV scale squark 
masses in models with minimal flavor violation. On the 
other hand, the neutralinos and charginos can naturally 
be light, near the weak-scale and with relatively small 
mass splitting. 

We consider production of a generic new heavy par- 
ticle Q in a collider, which decays in the following way: 
Q -> ■ ■ ■ X2 -> ■ ■ ■ XiX. Here xi,2 are the LSP and NLSP 
respectively, while X is a SM particle. The set of dots 
represents a multiplicity of SM states arising from the in- 
termediate stages of the decay chain. This rather general 
decay topology is shown in Fig. [T] If the mass difference 
between Q and \2 is large, mQ—m X2 ^> m X2 , the daugh- 
ter particle X2 will often be highly boosted. 

We point out that for a given X momentum configura- 
tion, there is a kinematic boundary for the missing mo- 
mentum. In the case of a boosted decay chain, the total 
missing momentum is constrained to lie in a narrow cone 
around the total X momentum, which we call the MET- 
cone. This observation is the central focus of this letter, 
and motivates the construction of a simple variable which 
contains information about the mass spectrum. We begin 
with a brief discussion of collinearity and its dependence 
on the mass parameters. 2 




FIG. 1: A schematic picture of a boosted decay chain. 

Collinearity and the MET- Cone — The kinemat- 
ics of two-body \2 decays are straightforward. We take 
the \2 particle to have relativistic boost factor 7 and 
velocity j3 in the lab frame. The angle X2 x between 
the visible particle X and the direction of motion of the 
parent \2 is then given by 



1 From here on, we use the semantics of supersymmetry, and re- 
fer to the lightest particle charged under the discrete symmetry 
as the LSP, and the next lightest as the NLSP. However, our 
method is generic, and applies to all TeV scale physics in which 
the exotica carry a conserved Z2 charge under which the SM 
fields are neutral. 



tan# X2 x = 
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7 V/¥ cos6» + /3 



(1) 



2 There is a recent study 1121 which also uses boosted decays for 
mass measurement, but in a quite different way. 
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where 9q is angle between X and \2 in the rest frame of 

X2, 



and /3q is the velocity of X in the rest frame of \2- 



The angle takes on values in the range (for /3q < 



< tanfl 



< 



X 



7>1 



7 



7/3 yTl - {pffpf 



(2) 



The velocity /3g is a function of the masses of the three 
particles involved, and it characterizes the allowed phase 
space of the X2 decay. The angle X 2Xi can be obtained 
by exchanging mx with m Xl in the above equations. A 
collinear configuration is achieved with a large 7, and 
with narrow phase space for the X2 decay. 

The \2 boost factor is determined by several variables. 
As a simple example, we consider a heavy exotic, with 
mass vtiq which decays to a massless SM particle (e.g. a 
jet) and the NLSP, with mass m X2 . For toq ~ 2 TeV 
and m X2 ~ 200 GeV, a boost factor of 7 = 5 is achieved 
in the rest frame of Q. However, at a hadron collider the 
Q particle will be produced with some transverse as well 
as longitudinal momentum, providing a distribution of 
boost factors. In addition, in multistage decay chains, the 
typical boost factors will depend on the mass spectrum 
of particles participating in the cascade. 

The boost factors of X and \i m the lab frame are 
given by 



7 Xl ,x=77o Xl ' X (l±/W' X costf ) 



(3) 



The magnitudes of the 3-momenta of X and xi m the 



Xl" l Xl 



and px = 



lab frame can be written as p Xl 
IxPxmx- 

Eq.'s (2j3l, define a kinematic boundary on the con- 
tribution of one Xi to the total fij>. These kinematic 
endpoints persist when there are two Xi particles in a 
single event. We illustrate this in Figure [2] where we 
display the region allowed for the total JSt vector in each 
event for given NLSP and LSP masses. We assume an 
event topology where all fHj> arises from two Xi particles 
(i.e. there are no neutrinos in the event), and that each 
of the two decay chains terminates with the NLSP decay- 
ing to a Z-boson plus the LSP. We restrict to a particular 
configuration of Z momenta. The j/-axis reflects the com- 
ponent of the ]£t vector parallel to the total transverse 
Z-momentum while the a;-axis displays the remaining ]£t 
vector component. As expected, the total ]£t vector is 
correlated with the Z-momenta, with kinematic bound- 
aries determined by the mass spectrum of the underlying 
theory. 

The m*£f variable — We utilize the collinear limit 
to inspire a test variable whose distribution yields the 
masses of the exotica. In the case of small mass splitting, 
Am = m X2 — m Xl — mx, the decay products are not 
significantly relativistic in the rest-frame of the parent, 
X2- Thus in the lab frame, where the X2 has relativistic 
velocity, the boost factors of all three particles are nearly 
the same, and the xi particles are closely aligned with 
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FIG. 2: In this figure we display the MET-cone boundary for 
different NLSP and LSP masses (keeping the mass splitting 
constant). The total $t vector must lie within the boundary 
for particular choices of exotica masses. The Z bosons are in 
a configuration where each has boost factor 5, both he in the 
transverse plane, and are separated by a 90 degree angle. 

the X-momenta. With this scenario as our motivation, 
we define "test" missing 3-momenta as 



Px 



,b m x? 
mx 



(4) 



with m^ 5 * defined for each event by minimization of the 
following quantity: 



, T.total 



24 



(5) 



This is an analytic procedure, as this formula is quadratic 



The minimization results are given by m 



test 
Xl 



m x ^T,y/PxT and A^ in /f T = fa J fx 



The vari- 



ables 



and A_^™ ln /_^T rescale respectively the y and 



x components of the J£t vector event-by-event. 

As a means of quality control, we veto signal events 
in which A|?™ ln / J£ T > e, for some sufficiently small e. 
The efficiency of such a cut is itself a rough measure of 
the mass splitting. In Figure [3j we show rescaled MET- 
cones in the m 1 ^ vs A|?™ m /|? T plane for both a scenario 
where the mass splitting is very small, and another where 
an O(10%) splitting is assumed. From these figures, one 



observes that 



Xi 



has two approximate endpoints in 



the small A^?™ n /^? T region, where the cones all inter- 



sect the 



\i 



axis. These endpoints are a manifestation 



of the MET-cone boundaries when projected onto the 
transverse plane, as illustrated in Figure [4] 

The limit e — > is equivalent to an alignment condition 
on the momenta: pJ} total — >. 0^ > total mt cst / m x- In this 



Xl 



limit, one can express m^ 5 ' in terms of the measurable 
parameters of the event. The result can be written as 
an expansion in the angular separation between X and 
Xi for both sides of the decay chain, 6 a,b , in the near- 
collinear case. In a configuration where both X's are in 
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FIG. 3: In these two plots, we compare the rescaled MET- 
cones of two scenarios with decay X2 — ► Xi^i with identi- 
cal X2 mass (200 GeV). In the top (bottom) plot, we take 
Am = m X2 - (m xi + m z ) = 1 GeV(10 GeV). Different 
contours correspond to randomly chosen Z-momentum con- 
figurations. Note that Tn!j, est has two endpoints in the small 
A|?™ m /|?T region. In addition, the very small mass splitting 
scenario (top) predominantly occupies a much narrower range 

Of Afi^/fiT. 



the transverse plane, we obtain a relatively clean result: 



m. 



7o Xl l 



cos 0g 



/3 cos ( 



x (1 - cot 6>f b cos </> a 6» a + esc 9* b cos 6 6 ) , (6) 

where /3 and 7 refer to the NLSP in the a-chain, and 
/3q ss /Jq 1 is used. Here (9 a , <f> a ) are the spherical coor- 
dinates of p^ in the lab frame where the z-axis is along 
the fx', 9ab 1S ^ e an gl e between p£ and p x . At zeroth 
order in the 9 a ' b expansion, the endpoints of m 
given by: 



tost 
Xl 



arc 



(7) 



which are achieved when 9q = and 7r respectively, i.e. 
Xi moving forward or backward along the \i boost di- 
rection in its rest frame. The range of m^ t becomes 
smaller as the phase space for the xi decay shrinks. In 
the relativistic limit, (3 ~ 1, the endpoint positions m± st 



approximately determine the unknown masses 
m Y , . 



and 




FIG. 4: Illustration of the MET cone and its relation to the 
mi?* variable. 



Now let us discuss the non-collinear corrections to 
the endpoint positions from Eq. First we note 

that 9 a depends on Oq and in the near-collinear case 
it does not shift the endpoint positions. The domi- 
nant contribution comes from 9 b which is independent 
of 9q and gives rise to a shift of the endpoint position 



Amf « ±m!£ st esc 9* b b . Since 9" follows a distribu- 
tion determined by the kinematics of the decay, it leads 
to a smearing of the distribution near these endpoints. 
In the near-collinear case, the variation of 9 b around the 
central value is small, and smearing is minimal. 

In the general case where the X's are not in the trans- 
verse plane, one needs to project the MET cone into the 
transverse plane and then impose the alignment condi- 
tion. One can still express the result in a collinear expan- 
sion. The zeroth-order result remains the same as that in 
Eq. Q. However, the higher order expansion coefficients 
are modified by trigonometric functions expected to be 
of order one. 

In summary, if the endpoint positions of the m^ 3 * dis- 
tribution are measured from the data, we can find so- 
lutions for the masses of the LSP and NLSP using the 
relation in Eq. Q. 

Numerical Results — We now explore the effec- 
tiveness of the m^ 5 ' method using Monte-Carlo simula- 
tion. We consider squark pair production qL<j_L followed 
by the decays ql — > qj(2 — > <IXiZ in SUSY models. We 
consider the four spectra shown in Table [I] For each of 



Model 


m xi 


m X2 


m<j L 


^^^test \ thco 


(mV cst ) theo 


1 


100 


200 


1000 


54.6 


183.2 


2 


100 


250 


1250 


21.6 


463.0 


3 


200 


300 


1000 


117.9 


339.2 


4 


200 


350 


1250 


52.6 


761.0 



TABLE I: The relevant masses in four SUSY models and the 



expected endpoints (m± 



test \ thco 



. Masses are given in GeV. 



these models, we simulate 20k events of squark pair pro- 
duction and decay in pp collisions at yfs = 14 TeV in 
MadGraph [T5] using the 2 — » 6 matrix element. Events 
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are selected according to the parton-level cuts in Table 
HH The cos 6 ab cut is to ensure the coefficients in the 9 a ' 
expansion are not too large, which would obscure the 
endpoints. For Models 2 and 4, we have slightly loos- 
ened the selection cuts in order to get better statistics 
near the tail of the distribution: a) e = 0.2 for Model 2 
and 4. b) |r^> tot | < 3.0 for Model 4. 















Pt 


\v z \ 


|^,tot| 


fir 


e 


C° S d ab\ 


> 50 GeV 


< 3.0 


< 1.0 


> 200 GeV 


0.15 


< 0.5 



TABLE II: Event selection cuts for m£f . Here rj z 
pseudorapidity of the total Z momentum. 



is the 




m<£ (GeV) 



FIG. 5: The distributions of m£f for Model 1 and 2. For 
Model 2, we have normalized the distribution by a factor of 
two. 

In Fig. 5j we show the m^ 5 * distribution for Models 1 
and 2. Models 3 and 4 are very similar. For Model 1, 
we can see that the distribution is approximately a trian- 
gle with two endpoints at around 50 and 200 GeV. The 
shape of the distribution near each endpoint is slightly 
smeared due to deviations from collinearity. A simple 
way to extract these endpoints can be achieved by a lin- 
ear fit and taking the x- intercept. This would typically 
under-estimate m!f st while over-estimating m+ st . In our 
analysis, we take the position at the half maximum near 
the lower edge for m^ st to get a better estimation. For 
the upper endpoint, we use the x-intercept of a linear 
fit. More complicated fits and estimation are certainly 
possible. 

The results are shown in Table |III| together with the 
statistical errors. The results are consistent with the ex- 
pected zeroth-order endpoints, as shown in Table [I] For 
Model 2, as seen in Fig. [5j the upper endpoint is much 
less sharp than the lower one. Estimation of its position 
is subject to a relatively large systematic uncertainty de- 
pending on the binning and the choice of the fitting re- 
gion. Fortunately, the calculated masses are not very 
sensitive to the upper endpoint position. For a reason- 
able estimate of the m+ st in the range 400 — 500 GeV, the 
calculated masses (m™ eos , m™ eas ) vary only mildly from 
(103 GeV, 241 GeV) to (116 GeV, 264 GeV). Therefore, 



Model 


m xl 


m X2 




mV cst 


mens 
" l Xl 


mens 
" l X2 


1 


100 


200 


55 ±2 


205 ±3 


106 ± 2 


208 ±3 


2 


100 


250 


27 ± 2 


454 ± 20 


110 ±5 


253 ±5 


3 


200 


300 


112 ±5 


342 ± 10 


195 ±5 


296 ±5 


4 


200 


350 


49 ±2 


682 ± 16 


183 ±5 


329 ±5 



TABLE III: Results of the measured m^ 1 endpoints and 
m xi 2 for four SUSY models. Masses are given in GeV. 

even in this relatively less collinear case we still obtain 
a good estimate of the masses. For all four models, the 
estimated endpoints and calculated masses are summa- 
rized in Table ITTT1 The fitted masses are all within ~ 10% 
of the true masses. 

Summary and Outlook — In this letter, we explored 
a novel method of measuring the absolute mass scale of 
exotic particles in events with missing energy which is 
inspired by boosted cascade decays. This method uses 
the fact that in the boosted decay there is a limited vari- 
ation in both the direction and magnitude of the total 
three-momentum of missing particles relative to the total 
three-momentum of the visible partners. The boundary 
of the allowed region, or the MET-cone, is determined by 
the mass parameters and the configuration of the visible 
particles. We constructed a variable m^ 5 ' which has has 
endpoints in the approximate collinear case. These end- 
points depend on the masses involved in the final step 
decay, and once observed from the data, can be used to 
determine these masses. 

The m^ 5 * variable works best in the collinear limit. 
Given the data, the evidence of collinearity in the final 
step decay can be seen in various ways. First, one would 
see a peak at in the A^?™ ln / fi T distribution. Second, 
one would see well-defined endpoints in the m^ 5 * distri- 
bution. Once the masses of xi an d Xi are measured, one 
can find the masses of heavier exotics upstream of the 
NLSP decay using more standard techniques. 

We have demonstrated our method in a setup with two 
symmetric decay chains with a two-step cascade decay, 
however it applies for longer decay chains as well. An 
advantage of this method is that one does not require 
information on all of the visible SM particles. 

The m^ 5 * variable has the advantage that it is simple 
to calculate on an event-by-event basis. However, it does 
not fully utilize the information available in the MET- 
cones. Developing a more effective method which takes 
full advantage of these kinematic boundaries is currently 
under investigation. 
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